Magnetic polarizability of hadrons from lattice QCD in the 

background field method 



Frank X. Lee and Leming Zhou 

Center for Nuclear Studies, Department of Physics, The George Washington University, 

Washington, DC 20052, USA 



Walter Wilcox 

Department of Physics, Baylor University, Waco, TX 76798, USA 

Joe Christensen 

Physics Department, McMurry University, Abilene, TX 79697, USA 



Abstract 

We present a calculation of hadron magnetic polarizability using the tech- 
niques of lattice QCD. This is carried out by introducing a uniform external 
magnetic field on the lattice and measuring the quadratic part of a hadron's 
mass shift. The calculation is performed on a 24 4 lattice with standard Wilson 
actions at beta=6.0 (spacing a = 0.1 fm) and pion mass down to about 500 
MeV. Results are obtained for 30 particles covering the entire baryon octet 
(n, p, S°, E", E+, S", 3°, A ) and decuplet (A , A", A+, A++, E*°, E*~, 
E*+, ~*°, 3*", Q-), plus selected mesons (vr°, vr+, vr~, K°, K+, K~, p°, p+, 
p~ , K*°, K* + , K*~). The results are compared with available values from 
experiments and other theoretical calculations. 

I. INTRODUCTION 



Electric and magnetic polarizabilities are fundamental properties of hadrons. They de- 
termine the dynamical response (deformation) of a hadron to external electromagnetic fields, 
and provide valuable and sensitive information about the internal strong interaction struc- 
ture of the hadron as a composite particle. In [1] , the results of a lattice calculation of hadron 
electric polarizability for neutral hadrons, based on the methods of [2], has been presented. 
In this work, we will focus on the first calculation of hadron magnetic polarizability from 
lattice QCD. Our aim is to present the lattice data in sufficient detail in order to facilitate 
comparison with existing information in this area. 

Traditionally, the symbol a is used to represent the electric dipole polarizability, and 
j3 for the magnetic dipole polarizability. A considerable number of experiments have been 
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performed with the aim of measuring nucleoli polarizabilities. Most of these have specif- 
ically targeted the proton, some attempts have also been made at measuring the neutron 
polarizabilities. The experimental values are listed in Tables I and II. 

TABLE I. A list of electric and magnetic polarizabilities in units of 10 -4 fm 3 for the proton 
from recent experiments. The method is elastic Compton scattering off the proton. 



Illinois (1991) [3] 10.9±2.2 ± 1.3 3.3=f2.2 =F 1.3 

Mainz (1992) [4] ^KiM+im ^i^im 

Saskatoon (1993) [5] 9.8±0.4 ±1.1 4.4^0.4^1.1 

Saskatoon (1995) [6] 12.5±0.8 ± 0.5 2.1+"0.8 T 0.5 

MAMI (2001) [7] 11.9±0.5=f 1.3 ± 0.3 1.2±0.7 ± 0.3 ± 0.4 



Despite some shortcomings, measurements for the proton polarizabilities in these exper- 
iments are in reasonable agreement with each other. The experimental situation regarding 
the polarizabilities of the neutron is still quite unsatisfactory. This is mainly because a 
neutron Compton scattering experiment cannot be directly performed. 

TABLE II. A list of electric and magnetic polarizabilities in units of 10~ 4 fm 3 for neutron from 
recent experiments. In the method column, Elastic means elastic (or free) Deuteron Compton 
Scattering; Quasifree means Quasifree Deuteron Compton Scattering. 



Ref. Method a n (3 n 

~\8] Neutron Scattering 12.3±1.5±2.0 3.1±1.6±2.0 

[9] Quasifree 12.5±1.8lJ;J ± 1.1 2.7=Fl.8l?;f T 1.1 

[10] Quasifree 7.6 ~ 14.0 1.2~7.6 

[11] Elastic 5.5±2.0 10.3^2.0 

[12] Elastic 8.8±2.4 ± 3.0 6.5 =f2.4 =f 3.0 



The experimental situation on the nucleon can be approximately summarized as the 
following: the electric polarizability is roughly the same for the proton and neutron, with a 
value of around 10 in units of 10 -4 fm 3 ; the magnetic polarizability is roughly the same for 
the proton and neutron, with a value of around 3 in the same units. 

On the theoretical side, nucleon polarizabilities have been most studied in the frame- 
work of chiral perturbation theory (ChPT) [13-18]. Other approaches include quark mod- 
els [19-21], chiral soliton models [22,23]. For reviews on polarizabilities, see Refs. [24,25]. 
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II. METHOD 



A. Weak field expansion 

The approach is based on the mass shifts of the particles measured both in the absence 
and presence of magnetic fields in the QCD vacuum. For small external magnetic fields, the 
mass shift (in the unit system of h — 1 = c) 

Am{B) = m{B) -m(0), (1) 

is given by 

Am(B) — —ft - B - -(3B 2 . (2) 
By averaging Am{B) over the field, B, and its inverse, —B, we will form 

Am(B) even = -\{3B 2 - (3) 

— * 

After we get the even B mass shift from the lattice simulation, we do least-chi-square fits 
to the data points with a polynomial 

Am(B) even = c 2 B 2 + c 4 B 4 + ■■■. (4) 

The magnetic polarizability is then the negative quadratic coefficient 

P = -2c 2 . (5) 

The quartic and higher terms in Eq. (4) are included in order to measure possible numerical 
contamination. 



B. Interpolating fields 

The calculation of hadron masses centers around the time-ordered, two-point correlation 
function in the QCD vacuum, projected to zero momentum: 

G(t)=£(0|T{^(0)}|0>, (6) 

X 

where 77 is the interpolating field built from quark fields with the quantum numbers of 
the hadron under consideration. On the quark level, G(t) is evaluated in terms of quark 
propagators by way of path integrals. On the hadronic level, G(t) is a sum of exponentials 
which has the entire mass spectrum of a given channel. The ground state can be extracted 
by fitting G(t) at large times where it dominates. 

We consider a wide variety of particles in this study. For octet baryons, the local inter- 
polating fields are 
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rf = e abc (d aT C-f 5 u b )d c , 



v p = t abc {u aT c lb d h )u c , 



(7) 
(8) 



rf =e abc (d aT Cj 5 s b )d c , 



so = l abc 
V2 



(u aT C lb s b )d c + (<f T C 75 s> 



aT/-< „b\ 



(9) 
(10) 



rf + = e abc (d aT C lb u b )s c , 



rf = e abc (s aT C-f 5 d b )s c , 



(11) 
(12) 



?f T = t abc {s aT C lb u b )s c , 



b\ n c 



(13) 



A 8 _ 1 abc 



2(u aT C l5 d b )s c + (u aT C l5 s b )d c - (d aT C l5 s b ) 



u 



(14) 



In the above expressions, C is the charge conjugation operator and the superscript T means 
transpose. The sum over the color indices is implied and the e abc ensures that constructed 
states are color- singlet. For simplicity, the explicit dependence of the quark field operators 
on space-time, q(x), is not written out. The normalization factors are chosen so that in the 
limit of SU(3)-flavor symmetry all correlation functions simplify to that of the proton. In 
addition to the octet lambda A 8 , we also consider flavor-singlet lambda A 5 , 



r/ AS = e abc e uds {u aT C l5 d b )s c . 
Using the transpose of the terms in the parenthesis, it may be written as 



r/ AS = -2t abc \-(u aT C-f 5 d b )s c + (u aT C-f 5 s b )d c - (d a T C 75 s t 



(15) 



(16) 



which has a structure similar to i] AS except for the coefficient of the first term. Since SU(3)- 
flavor symmetry is broken by the strange quark, it is interesting to study an interpolating 
field that is made up of the terms common to both types, 



AC = take 

V2 



[u aT C lb s b )d c - (d aT C^s b )u c 



(17) 



Note that since the u and d quarks respond to magnetic fields differently, the usual isospin 
symmetry in u and d quarks is explicitly broken. 

For decuplet baryons, the local interpolating fields are 



= e abc (rf aT C 7 ^ 6 )d c , 



(18) 



4 



A 



abc 



2(d aT C-f^u b )d c + (d aT C^d b )u ( 



(19) 



,A+ 



r abc 



2(u aT C^d b )u c + (u aT C^u b )d c 



(20) 



nf~ = -^e abc [2(d aT C^s b )d c + (d aT C^d b )s c 



(21) 
(22) 



[u aT C^d b )s c + {d aT C llx s b )u c + {s aT C ltl u b )d c 



(23) 



2(u aT C^s b )u c + {u aT C lfl u b )s c 



(24) 



a6c 



2( S aT CVV + (* aT CVV 



(25) 



r abc 



Hi 



(26) 



We use the Lorentz index \x = 3 in this calculation. Note that for baryons the correlation 
function G(t) is a 4x4 matrix in Dirac space. In the Dirichlet boundary condition used in 
this work, the sum of the upper diagonal components couple to the positive-parity state, 
while the lower components couple to the negative-parity state. So only diagonal elements 
are computed. 

For mesons, we consider the pseudo-scalar states 



rf ( x ) = u(x)~f 5 d(x), 



(27) 



if (x) = —= u(x)^ 5 u(x) + d(x)~f 5 d(x) 
v2 



(28) 



rf + (x) = d(x)~f 5 u(x), 



r] (x) = u(x)>y 5 s(x), 



r] K °(x) = d(x)-f 5 s(x), 



(29) 
(30) 
(31) 
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and the vector states 



r ] K+ (x) = s(x) l5 u(x), (32) 



% ( x ) = u(x)^d(x), (33) 



u(x)-i^u(x) + d(x)^^d(x) 



(34) 



ri? + (x) = d(x)^u(x) : 



% ( x ) = s(x)^d(x), 
V K *°( X ) = d( x h»s( x ), 



(35) 
(36) 
(37) 



V/j, ( x ) = s(x)7 M u(x). 



(38) 



Here q = q +r jo- For the vector mesons, we average over the spatial Lorentz indices fi = 1, 2, 3. 
Note that in forming the correlation function for the 7r° state, there is a term from the local 
contraction in the form of 



V2 



M- 1 (x,x)-Mj 1 (x,x) 



(39) 



where M~ 1 (x, x) denotes all-to-all quark propagators (disconnected quark loops). In the 
zero magnetic field and SU(2) isospin symmetry in u and d quarks, the term vanishes. In 
the presence of magnetic field, however, it is no longer so. The same is true for the p° state. 
The calculation of these disconnected loops are prohibitively expensive and is prone to large 
statistical errors. In this work we will ignore the effect of the disconnected loops for these 
mesons. 



C. Lattice Techniques 

We briefly discuss how the magnetic field is introduced on the lattice. The procedure we 
use is very similar to but not exactly the same as the one presented in [2]. 

In the continuum case, the fermion action is modified by the minimal coupling prescrip- 
tion 

d» -> <9 M + iqA^ (40) 

where q is the charge of the fermion field and is the vector potential describing the 
external field. On the lattice, the prescription amounts to supplying a phase factor, e taqAlJ -, 
to the link variables in a given direction. Choosing 
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A 2 = Bxi and A = A x = A 3 = 



(41) 



a constant magnetic field 

F 12 = d x A 2 - d 2 A 1 = B 3 = B (42) 

can be introduced in the z-direction. 

We introduce two dimensionless parameters which characterize the field. One is given 

by 

r] = qBa 2 . (43) 

The other is the integer lattice length, 

p = x 1 /a. (44) 

In terms of these two parameters the phase factor becomes: 

e iaqM = e ir,p ^ (1+iVP), (45) 

where we have linearized the phase factor to mimic the continuum prescription. The field 
strength we use will be chosen small enough to satisfy the linearization requirement [2]. 

To summarize, the method to place the external magnetic field on the lattice in the z- 
direction is to multiply each gauge field link variable in the y-direction with a re-dependent 
factor: 

U 2 (x) -> (l + ir) P )U 2 (x) (46) 



D. Simulation Details 

Most of our results are based on the standard Wilson quark action on a quenched 24 4 
lattice with lattice spacing a = 0.1 fm. The lattice coupling (3 = 6.0. We have analyzed 150 
configurations to extract hadron magnetic polarizabilities. Fermion propagators M -1 were 
constructed at six different quark masses, which correspond to six the k values, 

k = 0.1515, 0.1525, 0.1535, 0.1540, 0.1545, 0.1555. (47) 

The critical kappa value is k c = 0.157096. Using 

"* = Ta ( \ - b (48) 

we can get the corresponding quark masses are: 230 MeV, 189 MeV, 147 MeV, 126 MeV, 
105 MeV, 64 MeV. The corresponding pion masses are: 1000 MeV, 895 MeV, 782 MeV, 721 
MeV, 657 MeV, 512 MeV. The strange quark mass is set at k — 0.1535. The lattice source 
is (x,y,z,t)= (12, 1, 1, 2). 
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In an exploratory study, the same calculation was performed using a different action, 
namely, the tree-level tadpole-improved Liischer-Weisz gauge action with lattice spacing 
a = 0.17 fm (or 1/a — 1159 MeV) set from the string tension, and the tadpole-improved 
clover quark action. We accumulated 15 configurations and the preliminary results have been 
presented in Ref . [26] . We will focus on the Wilson results. The clover results will be shown 
together with the Wilson results mainly for comparison purposes and as an independent 
check. In most cases, the results from the two action schemes are consistent with each other 
with errors, giving us confidence in the calculation. 

We used six different values of the parameter in units of 1CT 3 

rj = 0.0, +0.36, -0.72, +1.44, -2.88, +5.76. (49) 

The rj values in this sequence are related by a factor of —2. Thus we are able to study the 
response of a hadron composed of both up and down (strange) quarks, whose charges are 
related by the same factor, to four different nonzero magnetic fields. In units of 10 _3 e _1 a -2 , 
the magnetic field takes the values -1.08, 2.16, -4.32, and 8.64. In physical units, the magni- 
tude of the weakest magnetic field is 2.46 xlO 13 tesla. This is a very strong magnetic field. 
On the other hand, in the sense of the mass shift, this is really not that strong. We can 
estimate the ratio of mass shift of proton to the mass of proton: 

^ = iM! = L6 o x 10" 5 . (50) 
m m 

Here we have used 2xl0~ 4 fm 3 as the value of f3 p . From this rough estimation we can see 
that the mass shift is very small even in such a strong magnetic field. 

III. RESULTS 

In this section, we present the numerical results of our lattice study. We calculated the 
magnetic polarizability for 30 charged and neutral hadrons from their quadratic response 
to the field. Since this is a first calculation of these quantities, we will show effective-mass 
plots for the computed mass shifts for all of the particles considered at one value (the 
weakest) of the magnetic field. They form the basis of our analysis and give an unbiased 
view of the quality of our data. The mass shifts are fitted to the polynomial form discussed 
in section II A from which the magnetic polarizability is extracted. We group the results 
according to the particle types: octet, decuplet, and mesons. 

A. Octet Baryons 

Fig. 1 shows the effective mass plot for the neutron and proton. There is reasonable 
plateau behavior for the neutron and our results are extracted from the time window of 12 
to 14. The plateau behavior for the proton is not as good as for the neutron. Fig. 2 shows 
the mass shifts for the neutron as a function of the magnetic fields. There is good parabolic 
behavior going through the origin, an indication that contamination from the linear term 

— * — * 

has been effectively eliminated by averaging results from B and —B. It is typical of all the 
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particles so we will not show more such plots. We tried a number of ways in fitting to the 
data. First, we did the fitting with the B 2 term only, or with the B 2 term plus the B A term. 
We found that the B A term is numerically small. This confirms that the magnetic fields we 
use are indeed weak. We also tried fits using only the two smallest fields in magnitude, or 
the three smallest, or all four values. They all give results consistent within errors bars. The 
results quoted below are mostly from fitting to the two smallest fields. 

For charged particles, there is the possibility of Landau levels on the order of \qB / (2M)\ 
in the presence of magnetic fields, where q and M are the charge and mass of the particle, 
respectively. It is a linear term that is not eliminated by the averaging procedure. Their 
effects only show up at very large times, larger than where we fit the data. We studied this 
issue by including a linear term in the fit and found that the results for C2 are essentially 
unaffected. 

Fig. 3 shows the extracted magnetic polarizability for the neutron and the proton as 
a function of pion mass squared 1 in physical units. Two sets of data are displayed, one 
from the Wilson action based on 150 configurations, one from the clover action based on 15 
configurations. The two sets are consistent within errors. The experimental value, which is 
roughly the same for neutron and proton with a large uncertainty, is also indicated. One 
can see that the proton results are consistent with the experimental value, while the neutron 
results, which have smaller errors than the proton, are much bigger. This difference between 
the proton and the neutron is one of the surprises of our calculation. We do not attempt a 
chiral extrapolation in this work because the theoretical ansatz from chiral effective theories 
is quite limited at the moment and the systematic uncertainties of the lattice results are 
not quantified. The reader can get an idea on the chiral behavior by looking at the quark 
mass dependence in the figures and tables. The observations we make here and below are 
based on the value at the smallest pion mass and the trend as a function of the quark mass, 
therefore should be taken as qualitative. For example, the proton magnetic polarizability 
has a weaker quark mass dependence than the neutron. At the chiral limit, the agreement 
with experiment for the proton is expected to remain, and the difference between the proton 
and the neutron is expected to grow. 

Fig. 4 shows the effective mass plot for the octet sigma states. There is reasonable plateau 
behavior and we fit in the range of 12 to 14. Fig. 5 shows the magnetic polarizability for the 
three octet sigma states. The S + results are similar to those for the proton and also suffer 
from large errors. The S° results are positive and the S~ are negative, both have much 
smaller errors compared to the S + . 

Fig. 6 shows the effective mass plot for the octet cascade states, and Fig. 7 shows the 
corresponding magnetic polarizability extracted from 12 to 14. Fig. 8 shows the effective 
mass plot for the octet lambda states, and Fig. 9 shows the corresponding magnetic polar- 
izability. In this case, only the Wilson results are available. One can see that the results 
for A 8 and A c are almost the same. The signal for flavor-singlet A s plateaus much earlier 
than the other octet states, and the signal is lost quickly. Our results are extracted from 
the range 5 to 7. 



The pion mass squared is proportional to the quark mass in QCD. 
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FIG. 1. Effective mass plots for the neutron (upper) and proton (lower) mass shifts in lattice 
units at the weakest magnetic field. The lines correspond to quark masses from the heaviest (circles) 
to the lightest (triangles). 
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B [e^a- 2 ] (1CT 3 ) B [e^ 2 ] (1(T 3 ) B [e^a" 2 ] (lCT 3 ) 

FIG. 2. Mass shifts of the neutron as a function of the magnetic field B at the six quark masses 

(heavy to light in the order of top left-to-right to bottom left-to-right). 
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FIG. 3. Magnetic polar izability of the neutron and proton as a function of in physical units. 
The solid symbols are the results from the Wilson action, while the empty symbols are from the 
clover action. The Wilson results are obtained from the time window of 12 to 14. The experimental 
value, which we take to be the same for neutron and proton, is indicated by the star. 
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FIG. 5. Magnetic polarizability of the octet sigma states in physical units. The solid symbols 
are the results from the Wilson action, while the empty symbols are from the clover action. The 
Wilson results are obtained from the time window of 12 to 14. 

The results of our calculation in the octet sector from the Wilson quark action are 
summarized in Table III. 



Judging from the values at the smallest pion mass, one can observe the following broad 
features in the octet sector. The positively-charged states (p and S~) have small and positive 
values around 2, but suffer from large errors. The charge-neutral states (n, S°, 5°, A 8 ) have 
similar values to each other around 12. The negatively-charged states (£~, H~) have similar 
values to each other around —16. The A s has a very weak quark mass dependence. The 
values are expected to change a little bit after chiral extrapolations. 

In the case of hyperons, quark model [20] gives 1.7 for S + and —1.7 for E~, compared 
with our values of about 2 and —16, respectively. ChPT at the 0(p 3 ) order [14] predicts 
a = 10/? and small and positive values across the octet, which is a different pattern from 
the results on the lattice. It would be interesting to see if the pattern receives significant 
corrections at the next order of ChPT. 

B. Decuplet Baryons 

Fig. 10 displays the effective mass plot for the Delta states. The plateaus form earlier 
for the decuplet states than the octet states, probably due to the fact that the decuplet 
states are heavier. Fig. 11 shows the magnetic polarizability for the four Delta states. The 
interesting result here is that the A ++ polarizability is large and negative (around -60 at the 
smallest pion mass), which means it is most easily deformed (decreasing mass) under the 
probing magnetic field. It is the largest magnetic polarizability we observed on the lattice. 
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FIG. 6. Effective mass plots for the octet cascade mass shifts in lattice units at the weakest 
magnetic field in the order of H° (upper) and H~ (lower). The lines correspond to quark masses 
from the heaviest (circles) to the lightest (triangles). 
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FIG. 7. Magnetic polar izability of the octet cascade states in physical units. The solid symbols 
are the results from the Wilson action, while the empty symbols are from the clover action. The 
Wilson results are obtained from the time window of 12 to 14. 
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FIG. 8. Effective mass plots for the lambda mass shifts in lattice units at the weakest magnetic 
field in the order of A 8 (top), (middle), A s (bottom). The lines correspond to quark masses 
from the heaviest (circles) to the lightest (triangles). 
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FIG. 9. Magnetic polarizability of the octet lambda states in physical units from the Wilson 
action. The results for A 8 and hP are obtained from the time window of 12 to 14, while for hr 5 



to 7. 



TABLE III. The calculated magnetic polarizabilities for the octet baryons as a function of the 
pion mass from the Wilson action. The pion mass is in GeV and the magnetic polarizability is in 
10 -4 fm 3 . The time window from which each polarizability is extracted is given in the last column. 
The errors are statistical. 



K 


0.1515 


0.1525 


0.1535 


0.1540 


0.1545 


0.1555 


fit 


range 




1.000 


0.895 


0.782 


0.721 


0.657 


0.512 






p 


0.09 ± 0.29 


0.14 ± 0.37 


0.26 ± 0.48 


0.40 ± 0.56 


0.64 ± 0.67 


2.36 ± 1.20 


12- 


14 


n 


9.4 ± 0.2 


10.4 ± 0.3 


11.6 ± 0.4 


12.3 ± 0.5 


13.4 ± 0.6 


17.0 ± 1.1 


12- 


14 


£+ 


-0.15 ± 0.36 


0.09 ± 0.42 


0.24 ± 0.50 


0.40 ± 0.56 


0.61 ± 0.64 


1.60 ± 1.00 


12- 


14 


S° 


8.0 ± 0.3 


8.5 ± 0.3 


9.1 ± 0.4 


9.6 ± 0.5 


10.1 ± 0.6 


11.9 ± 0.9 


12- 


14 


£- 


-11.8 ± 0.3 


-12.6 ± 0.3 


-13.6 ± 0.4 


-14.2 ± 0.4 


-14.7 ± 0.4 


-16.1 ± 0.5 


12- 


14 




10.7 ± 0.3 


11.3 ± 0.4 


11.9 ± 0.4 


12.3 ± 0.4 


12.8 ± 0.5 


13.9 ± 0.7 


12- 


14 




-12.6 ± 0.3 


-13.1 ± 0.3 


-13.8 ± 0.4 


-14.1 ± 0.4 


-14.6 ± 0.4 


-15.6 ± 0.5 


12- 


14 


A 8 


9.1 ± 0.4 


9.9 ± 0.5 


10.8 ± 0.6 


11.4 ± 0.7 


12.1 ± 0.8 


14.0 ± 1.2 


12- 


14 


A c 


9.3 ± 0.4 


10.2 ± 0.5 


11.2 ± 0.6 


11.9 ± 0.7 


12.6 ± 0.9 


15.0 ± 1.4 


12- 


14 


A s 


3.6 ± 0.1 


3.4 ± 0.2 


3.3 ± 0.2 


3.3 ± 0.3 


3.2 ± 0.4 


3.2 ± 1.0 


5-7 
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FIG. 10. Effective mass plots for the Delta states mass shifts in lattice units at the weakest 
magnetic field in the order of, from top down, A ++ , A + , A , A - . The lines correspond to quark 
masses from the heaviest (circles) to the lightest (triangles). 
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FIG. 11. Magnetic polar izability of the A's in physical units. The solid symbols are the results 
from the Wilson action, while the empty symbols are from the clover action. The Wilson results 
are obtained from the time window of 10 to 12. 



Fig. 12 displays the effective mass plot for the decuplet sigma states, and Fig. 13 shows 
the corresponding magnetic polarizability as a function of the pion mass squared. Fig. 14 
displays the effective mass plot for the decuplet cascade states, and Fig. 15 shows the 
corresponding magnetic polarizability as a function of the pion mass squared. 

The results of our calculation in the decuplet sector from the Wilson action are summa- 
rized in Table IV. 

From the values at the smallest pion mass, one can observe the following features in the 
decuplet sector. The values are expected to change a little bit after chiral extrapolations. 
The positively-charged states (A + and £* + ) have negative values around -5 and relatively 
small errors. The situation is opposite to that for the octet members (p and The 
charge-neutral states (A , S*°, have similar values to each other around 10, similar to 
the situation in the octet sector. The negatively- charged states (E*~, Q~) have similar 
values to each other around -15, again similar to the situation in the octet sector. Note 
that the value for Q~ , -12.4(2), is simply predicted from the A" at the strange quark mass 
(k=0.1535). This prediction is free of uncertainty from chiral extrapolations. An experi- 
mental measurement of the Q~ magnetic polarizability is greatly desired. Historically, the 
precise determination of the mass and magnetic moment of the Q~ has paid significant div- 
idends to our understanding of hadron structure and dynamics. Its magnetic polarizability 
provides another such potential opportunity. 

There is no experimental measurements on the decuplet states, and very limited theo- 
retical information. An estimate in the Skyrme model [23] gives a value of -14.9 for A - , 
which is consistent with the value on the lattice. The large and negative value for the A ++ 
calls for theoretical explanations. 
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FIG. 12. Effective mass plots for the decuplet sigma mass shifts in lattice units at the weakest 
magnetic field in the order of E* + (top), S*° (middle), S*~ (bottom). The lines correspond to 
quark masses from the heaviest (circles) to the lightest (triangles). 
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FIG. 13. Magnetic polar izability of the decuplet sigma states in physical units. The solid 
symbols are the results from the Wilson action, while the empty symbols are from the clover 
action. The Wilson results are obtained from the time window of 10 to 12. 



TABLE IV. The calculated magnetic polarizabilities for the decuplet baryons as a function of 
the pion mass from the Wilson action. The pion mass is in GeV and the magnetic polarizability 
is in 10~ 4 fm 3 . The time window from which each polarizability is extracted is given in the last 
column. The errors are statistical. 



K 


0.1515 


0.1525 


0.1535 


0.1540 


0.1545 


0.1555 


fit 


range 




1.000 


0.895 


0.782 


0.721 


0.657 


0.512 






A++ 


-39.9 ± 0.7 


-43.9 ± 0.8 


-48.7 ± 1.0 


-51.5 ± 1.1 


-54.7 ± 1.3 


-63.1 ± 1.9 


10- 


12 


A+ 


-2.5 ± 0.2 


-3.0 ± 0.3 


-3.5 ± 0.5 


-3.9 ± 0.6 


-4.3 ± 0.7 


-5.1 ± 1.1 


10- 


12 


A 


7.6 ± 0.2 


8.2 ± 0.3 


8.8 ± 0.4 


9.2 ± 0.5 


9.6 ± 0.6 


10. 9± 1.0 


10- 


12 


A~ 


-10.1 ± 0.2 


-11.1 ± 0.2 


-12.4 ± 0.2 


-13.1 ± 0.3 


-14.0 ± 0.3 


-16.2 ± 0.5 


10- 


12 


£*+ 


-2.9 ± 0.3 


-3.2 ± 0.4 


-3.6 ± 0.5 


-3.9 ± 0.6 


-4.2 ± 0.6 


-5.1 ± 0.9 


10- 


12 


S*° 


7.9 ± 0.2 


8.4 ± 0.3 


8.9 ± 0.4 


9.2 ± 0.5 


9.5 ± 0.6 


10.3 ± 0.8 


10- 


12 


£*- 


-10.9 ± 0.2 


-11.7 ± 0.2 


-12.6 ± 0.3 


-13.1 ± 0.3 


-13.7 ± 0.3 


-15.0 ± 0.4 


10- 


12 


^*o 


8.1 ± 0.3 


8.5 ± 0.3 


9.0 ± 0.4 


9.2 ± 0.5 


9.4 ± 0.5 


9.8 ± 0.7 


10- 


12 




-11.9 ± 0.2 


-12.4 ± 0.2 


-12.8 ± 0.3 


-13.1 ± 0.3 


-13.4 ± 0.3 


-14.0 ± 0.3 
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FIG. 14. Effective mass plots for the decuplet H*° (upper), H*~ (lower) mass shifts in lattice 
units at the weakest magnetic field. The lines correspond to quark masses from the heaviest (circles) 
to the lightest (triangles). 
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FIG. 15. Magnetic polarizability of the decuplet cascade states in physical units. The solid 
symbols are the results from the Wilson action, while the empty symbols are from the clover 
action. The Wilson results are obtained from the time window of 10 to 12. 
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C. Mesons 



Fig. 16 displays the effective mass plot for the pion states. The plateau happens much 
later than in the neutron. We fit in the range of 16 to 18. In our calculation, n + and ir~ 
have identical mass shifts, meaning they have identical magnetic polarizability. The same 
is true for all other charged mesons we considered. This is expected because they are anti- 
particles to each other and have the same mass. Fig. 17 shows the corresponding magnetic 
polarizability as a function of the pion mass squared. 

Fig. 18 displays the effective mass plot for the kaon states. The quality of the plateaus 
is about the same as that of the pion states and we fit in the same range of 16 to 18. Fig. 19 
the corresponding magnetic polarizability as a function of the pion mass squared. 

Fig. 20 displays the effective mass plot for the rho mesons. The plateau for p° is not 
as good as that for it , and we fit in the range of 9 to 11. We also fit a bit earlier (range 
14 to 16) for p ± than for the pions (16 to 18). Fig. 21 shows the corresponding magnetic 
polarizability as a function of the pion mass squared. 

Fig. 22 displays the effective mass plot for the rho mesons, and Fig. 23 shows the corre- 
sponding magnetic polarizability as a function of the pion mass squared. 

The results of our calculation in the meson sector from the Wilson action are summarized 
in Table V. 

Based on the values at the smallest pion mass, one can observe the following features in 
the meson sector. The charged states (7r ± , K^) have similar values to each other around 
—24, which are about twice in magnitude (around —12) as the other two pairs (p^, K*^). 
The neutral states (p°, K°, K*°) have similar values to each other around +5, except n° 
which is about 3 times as large (around 15). 

The most recent measurement of pion polarizability [27] gives (a — (3) n + = 11.6, which 
implies /3 W + of about —6 assuming (a+/3) n + pa according to leading-order ChPT. The lattice 
value of —24 has the same sign but is much more negative. There is no measurement for 
the rho and K mesons, although there are plans to do so [28]. Other theoretical calculations 
have produced various values, about —2 for 7r ± , and a small but positive value (about +2) 
for 7T° They include ChPT at two-loop order [29], dispersion sum rule (DSR) [30], and NJL 
model [31]. The positive sign for tt° is consistent with the lattice value. 

In the case of kaon, ChPT predicts —0.6 for K + [32]. NJL model predicts —11 for K + 
and 13 for K° [33]. In the case of vector mesons, the information is extremely limited. The 
only calculation is from the QCD string theory [34] which gives —0.8 for rho and —0.6 for 
K*. 



IV. SUMMARY AND OUTLOOK 

In this work we have presented the results of the first lattice calculation of hadron 
magnetic polarizabilities. Hadron masses were extracted from fits to baryon and meson two- 
point correlation functions at six different quark masses and four different nonzero magnetic 
fields. Magnetic polarizabilities of hadrons were extracted from the hadron mass shifts 
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FIG. 16. Effective mass plots for the pion mass shifts in lattice units at the weakest magnetic 
field in the order of 7r° (upper) and tt^ (lower). The lines correspond to quark masses from the 
heaviest (circles) to the lightest (triangles). 
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FIG. 17. Magnetic polarizability for the pion states in physical units. The solid symbols are the 
results from the Wilson action, while the empty symbols are from the clover action. The Wilson 
results are obtained from the time window of 16 to 18. 
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FIG. 18. Effective mass plots for the kaon mass shifts in lattice units at the weakest magnetic 
field in the order of K° (upper) and (lower). The lines correspond to quark masses from the 
heaviest (circles) to the lightest (triangles). 
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FIG. 19. Magnetic polar izability for the kaon states in lattice units. The solid symbols are the 
results from the Wilson action, while the empty symbols are from the clover action. The Wilson 
results are obtained from the time window of 16 to 18. 
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FIG. 20. Effective mass plots for the rho meson mass shifts in lattice units at the weakest 
magnetic field in the order of p° (upper) and /r* 1 (lower). The lines correspond to quark masses 
from the heaviest (circles) to the lightest (triangles). 
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FIG. 21. Magnetic polarizability for the rho meson states in lattice units. The solid symbols 
are the results from the Wilson action, while the empty symbols are from the clover action. The 
Wilson results are obtained from the time window of 9 to 11 for p° and 14 to 16 for p^. 
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FIG. 22. Effective mass plots for the K* mass shifts in lattice units at the weakest magnetic 
field in the order of K*° (upper) and K*^ (lower). The lines correspond to quark masses from the 
heaviest (circles) to the lightest (triangles). 
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FIG. 23. Magnetic polarizability for the K* states in physical units. The solid symbols are the 
results from the Wilson action, while the empty symbols are from the clover action. The Wilson 
results are obtained from the time window of 14 to 16. 
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TABLE V. The calculated magnetic polarizabilities for the selected mesons as a function of 
the pion mass from the Wilson action. The pion mass is in GeV and the magnetic polarizability 
is in 10~ 4 fm 3 . The time window from which each polarizability is extracted is given in the last 
column. The errors are statistical. 



K 


0.1515 


0.1525 


0.1535 


0.1540 


0.1545 


0.1555 


fit range 


77V 


1.000 


0.895 


0.782 


0.721 


0.657 


0.512 




7T± 


-16.0 ± 0.3 


-17.8 ± 0.4 


-19.8 ± 0.5 


-21.1 ± 0.5 


-22.5 ± 0.6 


-26.4 ± 0.8 


16-18 


7T° 


8.4 ± 0.2 


9.1 ± 0.3 


10.1 ± 0.4 


10.8 ± 0.5 


11.6 ± 0.5 


14.1 ± 0.8 


16-18 


K± 


-18.4 ± 0.4 


-19.4 ± 0.4 


-20.5 ± 0.5 


-21.1 ± 0.5 


-21.7 ± 0.6 


-23.1 ± 0.6 


16-18 


K° 


3.7 ± 0.1 


3.8 ± 0.2 


4.0 ± 0.2 


4.1 ± 0.2 


4.3 ± 0.2 


4.7 ± 0.3 


16-18 


P ± 


-11.8 ± 0.3 


-12.5 ± 0.3 


-13.1 ± 0.4 


-13.3 ± 0.5 


-13.3 ± 0.6 


-12.5 ± 1.2 


14-16 


P° 


5.3 ± 0.2 


5.6 ± 0.2 


5.9 ± 0.3 


6.0 ± 0.4 


6.1 ± 0.4 


6.5 ± 0.7 


9-11 


K *± 


-13.0 ± 0.3 


-13.2 ± 0.4 


-13.3 ± 0.5 


-13.3 ± 0.5 


-13.2 ± 0.6 


-12.7 ± 0.8 


14-16 


K*° 


3.0 ± 0.2 


3.2 ± 0.2 


3.3 ± 0.3 


3.5 ± 0.3 


3.6 ± 0.3 


4.1 ± 0.5 


14-16 



induced by the presence of an external uniform magnetic field. The bulk of our results are 
based on 150 configurations on a 24 4 lattice at a spacing of about 0.1 fm and pion mass 
down to about 500 MeV, using the standard Wilson actions. The clover quark action was 
also used as a check. In both cases, clear signals for hadron magnetic polarizabilities are 
seen and the calculated results from Wilson and clover quark actions are generally in good 
agreement with each other. We investigated 30 particles sweeping through the baryon octet 
and decuplet and selected mesons. Most of the polarizabilities have not been measured, 
except for the nucleon and pion, so most of our results are predictions. Aside from the 
figures, the results are tabulated in three separate tables. Below is a summary of the main 
results of the calculation. 

Our value for the proton magnetic polarizability agrees reasonably with the most recent 
world average value of about 3 in units of 10~ 4 fm 3 , but suffers from relatively large errors. 
Our value for the neutron magnetic polarizability (about 15 to 20), which has smaller errors, 
is much greater than the world average value. This large difference between the proton and 
the neutron on the lattice is an interesting result that is worth further study. 

In the decuplet sector, the most interesting result is the large and negative magnetic 
polarizability (about —60) for the A ++ . Theoretical explanations, such as those from chiral 
effective theories, are called for. The tt~ value of —12.4(2) is a prediction free of uncertainty 
from chiral extrapolations, and can be directly compared with experiment. A measurement 
of the Q~ magnetic polarizability is greatly desired. 

Between the octet and the decuplet, we observe the following pattern. Positively-charged 
p and S + have relatively small and positive values, albeit with large errors, a situation 
opposite to that for the decuplet members A + and S* + which have negative and slightly 
larger values along with smaller errors. It is important to increase the statistics to obtain 
a better signal for the p and S + on the lattice to confirm this difference. In addition, the 
charge- neutral particles have values on the order of 20, while the negatively-charged particles 
on the order of —20. 

In the meson sector, we confirmed the expected result that a positively-charged particle 
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has identical magnetic polarizability as its negatively-charged partner (n"* 1 , K ± , p ± 1 K^). 
All charged mesons have negative magnetic polarizabilities, and all neutral mesons have 
positive ones. In terms of magnitude, the pseudoscalar mesons have about twice as large 
values as the vector mesons. Futhermore, the vector mesons have a weaker quark mass 
dependence than the pseudoscalar mesons. 

Taken as a whole, our results demonstrate the efficacy of the methods used in comput- 
ing the magnetic polarizabilities on the lattice. In addition to increasing statistics to the 
300 to 500 configurations range, further studies should focus on assessing the systematic 
uncertainties. 

First, finite-volume effects must be investigated by changing the lattice size or lattice 
spacing. This is particularly relevant considering the fact that the dipole magnetic polar- 
izability is a volume effect (in units of fm 3 ). Calculations on a larger lattice 32 4 are under 
way. Another volume is needed to perform a continuum extrapolation. 

Next, a chiral extrapolation is needed to provide better comparison with experiment and 
other approaches. To facilitate the chiral extrapolation, on the one hand, we need to perform 
calculations at smaller pion masses. This is hard to achieve with the standard Wilson quark 
action due to technical difficulties. An attractive alternate is the so-called twisted-mass 
quark action [35] which can push the pion mass down to about 250 MeV at a cost of about 
factor of two. To further decrease the pion mass, one needs chiral fermions such as the 
overlap [36] which has been used to push the pion mass down to about 180 MeV [37]. The 
downside of overlap fermions is its cost: it is much more expensive to simulate than the 
standard Wilson fermions. On the other hand, we need better ansatz from chiral effective 
theories to provide a physical basis for the extrapolations. 

Third, the effects of the quenched approximation should be investigated. Work is 
currently under way using the CP-PACS dynamical configurations [38] based on the 
renormalizat ion-group improved gauge action and mean-field improved clover quark action. 
The three lattice sizes in the set are useful for doing a continuum extrapolation. 

Finally, as far as the cost of our calculation is concerned, it is equivalent to 11 standard 
mass-spectrum calculations using the same action (5 values of the parameter 77 to provide 
4 non-zero magnetic fields, 5 to reverse the field, plus the zero-field to set the baseline). 
This factor can be reduced to 7 if only two non-zero values of magnetic field are desired. 
Reversing the field is well worth the cost: the magnetic moments can be extracted from the 
linear response in the mass shifts in the same data set [39]. 
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